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HR8799 is a A Bootis, 7 Doradus star hosting a planetary system and a debris disk with two rings. This makes this 
system a very interesting target for asteroseismic studies. In particular, this work is devoted to the determination of the 
internal metallicity of this star, linked with its A Bootis nature, and its age, taking the advantage of its 7 Doradus-type 
pulsations. To do so we have used the equilibrium code CESAM and the non-adiabatic pulsational code GraCo. We have 
applied the Frequency Ratio Method and the Time Dependent Convection theory to estimate the mode identification, 
the Brunt- Vai'sala frequency integral and the mode instability, making a selection of the possible models fulfilling all 
observational constraints. Using the position of the star in the HR diagram, the solar metallicity models is discarded. 
This result contradicts one of the main assumptions of the most accepted hypothesis explaining the A Bootis nature, the 
accretion/diffusion of gas from a star with solar metallicity. Therefore, in sight of these new results, a revision of this 
hypothesis is needed. The inclusion of accurate internal chemical mixing is necessary. The use of the asteroseismological 
constraints provides a very accurate determination of the physical characteristics of HR 8799: an age in the ranges [1 123, 
1625] and [26, 430] Myr, and a mass in the ranges [1.32, 1.33] and [1.44, 1.45] M , respectively, depending on the visual 
angle i. The determination of this angle and more accurate multicolor photometric observations can definitively fix the 
mass, metallicity and age of this star. In particular, an accurate age estimation is needed for a correct understanding of the 
planetary system. In our study we have found that the age widely used for modelling the system is unlike. 
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IT) ■ 1 Introduction 



Accurate ages stars hosting extrasolar planets are useful for, 
e.g., obtaining theoretical constraints on tidal interactions 
between Hot Jupiters and their hosts, dynamical modelling 
of exoplanetary systems, predicting exoplanetary luminosi- 
ties, etc, since they constrain the age of the planets. The 
eventual aim of refining the age determination of the stars is 
to obtain more information on planetary system formation 
to discern between the different scenarios. Several meth- 
ods can be applied to determine ages of stars, like the use 
of chromospheric activity proxies, lithium depletion, evolu- 
tionary track fitting, etc (Barrado y Navascues 1998; Bar- 
rado y Navascues et al. 1999). Nonetheless, for pulsating 
stars, asteroseismology is the best and most accurate method| 
we know that provides mass, radii and ages with an accu- 
racy order of magnitudes larger than obtained with other 
methods in the literature (Stello et al. 2009). 

The discovery by Marois et al. (2008), of the first plane- 
tary system by direct imaging around HR 8799 has made the 
age, mass and metallicity determination of this star a very 
important task. That was the starting point of a set of ap- 
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proximately ten studies about this system during 2009, none 
of them from an asteroseismic point of view, even consid- 
ering that there are three detected pulsational frequencies 
(Zerbi et al. 1999) that can be used to better understand the 
host star. 

The main characteristics of this system are: 

- The host star is a A5 Main Sequence star (Gray & Kaye 
1999), 7 Doradus pulsator (Zerbi et al. 1999), and it has 
A Bootis-type surface chemical peculiarities (Sadakane 
2006). 

- There are three sub-stellar objects orbiting the host star. 
They have been observed using direct imaging (Marois 
et al. 2008). 

- The system has a debris disk with two rings, one inner 
and another further out the three orbiting objects (Chen 
et al. 2009). 

- From a dynamical point of view, the stability of the sys- 
tem is a challenge, regarding the mass given for the ob- 
jects and star, and the age assumed for the system (Rei- 
demeister et al. 2009). 

These four characteristics make the planetary system 
HR8799 a very interesting one, with many challenges to 
understand its structure, formation and evolution. In partic- 
ular, the age of the star (and consequently of the system) 
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is one of the key quantities for all the studies done around 
HR8799, since it is directly related with the mass of the 
orbiting objects, and it imposes a goal for the dynamical 
stability of the system. 

In two recent works (Moya et al., 2010a, b), a compre- 
hensive modelling of HR8799 using asteroseismology has 
been done. The main conclusions are the following: 

- The use of asteroseismological observational data im- 
poses large constraints to the acceptable models of stel- 
lar internal structure. 

- The actual uncertainty in the knowledge of the age of 
HR 8799 is larger than shown in the literature. Only less 
than 20% of the accepted model using all the observa- 
tional constraints have an age in the range estimated by 
other works. 

- The internal subsolar metallicity obtained for the star 
provides some keys to understand the A Bootis nature. 

These works have shown the enormous potential of the 
use of asteroseismology to understand this system and oth- 
ers. Nevertheless, these studies cannot be conclusive due to 
the poor asteroseismic observational available data. There 
is not very accurate multicolour photometric observations. 
Therefore, additional observations are imperative to provide 
better asteroseismological constraints. 

2 Observational data 

The A5 V star HR 8799 has been extensively studied in the 
last years. The first studies of this star were done in the con- 
text of asteroseismology. Schuster & Nissen (1986) firstly 
reported HR 8799 as a possible SX Phoenicis type. Zerbi et 
al. (1999) observed this star in a multisite multicolour pho- 
tometric campaign, with Stromgren filters, and found three 
independent pulsational frequencies (/i = 1.9791 cd _1 , 
h = 1-7268 cd" 1 , and / 3 = 1.6498 cd" 1 , units are 
cycles per day), making it one of the 12 first 7 Doradus pul- 
sators known (Kaye et al 1999). This pulsating stellar group 
is composed of Main Sequence (MS) stars in the lower part 
of the classical instability strip. Their pulsation modes have 
periods in the range [0.5,3] days, that is, they are asymptotic 
g-mode pulsators. Gray & Kaye (1999) obtained an optical 
spectrum of HR 8799, and assigned an spectral type of kA5 
hFO mA5 V A Bootis, reporting an atmospheric metallic- 
ity of [M/H]=— 0.47 and accurate values of the stellar lumi- 
nosity, T e ff, and log g (see Table 1) . They also noted that 
HR 8799 may be also a Vega-type star, characterized by a far 
IR excess due to a debris disk. Sadakane (2006) developed 
a deep study of the metal abundances of this star, confirm- 
ing its A Bootis nature (with surface chemical peculiarities, 
Paunzen 2003). HR8799 is one of the three A Bootis stars 
with 7 Doradus pulsations (Rodriguez et al. 2006a, b) 

The use of different bolometric corrections changes sig- 
nificantly the value of the luminosity of the star. Using the 
Virtual observatory tool VOSA (Bayo et al. 2008), models 



Table 1 Physical characteristics of HR 8799 



T eff (K) 


7430±75 




logs (cm s -2 ) 


4.35±0.05 


Gray & Kaye 1999 


L(L ) 


4.92±0.41 




v sin i (km s _i ) 


37.5±2 


Kaye & Strassmeier 1998 


7r (mas) 


25.38±0.85 


van Leeuwen 2007 



with realistic metallicities best fitting the observations pro- 
vide a luminosity only 0.1 L/Lq larger than that given in 
Table 1 and, therefore, within the errors. On the other hand, 
the use of different parallaxes in the literature changes the 
value of the absolute magnitude of this star by less than 0.1 
L/Lq (see Moya et al. 2010a for details). 

3 Modelling of the observational data 

We have developed a grid of equilibrium models obtained 
with the CESAM code (Morel & Lebreton 2008). We vary 
the mass (in the range [1.25, 2.10] M with steps of 0.01 
M ), the metallicity (with values [M/H]=0.08, -0. 12, -0.32,| 
and —0.52), the Mixing-Length parameter MLT (values 0.5, 
1, and 1.5), and the overshooting (values 0.1, 0.2, and 0.3). 
The internal metallicity has been regarded as a free param- 
eter due to the A Bootis nature of the star, which hides its 
internal abundances. The mass, estimated to be 1.47 ± 0.3 
Mq by Gray & Kaye 1999, has been also regarded as a free 
parameter since it has not been directly determined. 

The pulsational analysis has been developed using the 
non-adiabatic pulsational code GraCo (Moya et al. 2004; 
Moya & Garrido 2008). Both codes have been tested in the 
work of the ESTA activity (Evolution and Asteroseismic 
Tools Activities, Lebreton et al. 2008; Moya et al. 2008). 
Using these tools, we performed a theoretical study of HR 8799| 
in an attempt at constraining physical and theoretical param- 
eters. In this work we will follow the same scheme used for 
the study of the 5 Scuti pulsators RV Arietis, 29 Cygnis and 
9 Aurigae (Casas et al. 2006; Casas et al. 2009; Moya et al. 
2006). Taking advantage of the 7 Doradus pulsations, we 
have used the Frequency Ratio Method (Moya et al. 2005) 
to estimate the possible mode identification of the observed 
modes, and the Time Dependent Convection (Gricahcene et 
al. 2005) for the analysis of the mode instability and the pos- 
sible spherical order comparing with observed multicolour 
photometry. All these techniques provide additional con- 
straints to the position of the star in the HR diagram, re- 
ducing the possible physical characteristics of the models 
fulfilling all the observations. 

4 On the A Bootis nature of HR 8799 

The first selection of models was done imposing only the 
spectroscopic constraints displayed in Table 1 . This selec- 
tion shows that there are no models with solar metallicity 
fulfilling observations, since the stellar luminosity derived 
is smaller than any of the possible luminosities of models 
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Fig. 1 T ff - Age diagram for the models of our grid ful- 
filling the spectroscopic observations, i.e., T e //, logg and 
luminosity displayed in Table 2. The age range obtained is 
[10,2337]Myr. 

with solar metallicity. This contradicts the main assumption 
of the theories explaining the A Bootis nature, i.e. the accre- 
tion/diffusion scenario where these stars have solar metal- 
licity, whereas the observed abundances are due to surface 
phenomena (Turcotte & Charbonneau 1993). 

The result of our work, discarding solar metallicity as 
the internal metallicity of the star, together with the fact that 
some A Bootis stars have debris disks not connected with the 
star (Chen et al. 2009), makes this accretion/diffusion sce- 
nario unlikely, but not negligible (Su et al. 2009). Therefore, 
the study of internal chemical mixing processes seems to be 
the key to explain the A Bootis nature, at least for HR 8799, 
as the solar abundances for C, N, O and S observed on its 
surface have still to be explained. 

5 On the mass and age of HR 8799 

There are several age determinations of the age of HR 8799 
in the literature (see Moya et al. 2010b for a complete re- 
view), most of them estimating its age in a range around 
[30,160] Myr. The most roboust determination found is the 
kinematic age determination (UVW). But the authors warned| 
that this method is not always reliable. 

The main complementary argument to that statistical de- 
termination is the position of the star in the HR diagram. 
In Moya et al. 2010b we have demostrated that this proce- 
dure does not provide accurate age estimations due to the A 
Bootis nature of HR 8799 (Fig. 1). This result was previusly 
pointed out by Song et al. (2001), how stimated the age of 
HR8799 in the range [50,1128] Myr. In Section 4 we saw 
that this nature hides the real internal metallicity of the star. 
The main consequence of this result is that the models ful- 
filling observations are in a range of ages [10,2337] Myr, 
a much broader range that one estimated by other authors 
of [30,160] Myr (Marois et al. 2008). Only a small amount 
(18.1%) of models in our representative grid have ages in 
the range mainly claimed in the literature. 



Fig. 2 T c ff - Age diagram for the models of our grid ful- 
filling the spectroscopic observations plus all the astero- 
seismological constraints (FRM + multicolour photometry 
+ instability analysis, see Moya et al. 2010a, b) for a rota- 
tion velocity V Iot ~ 60km s _1 . The age ranges obtained are 
[1123,1625], [26,430] Myr 

Therefore we need aditional constraints for an accurate 
estimation of the age and mass of this star. We have ex- 
plained in Section 3 that one of the techniques used for 
this study is the Frequency Ratio Method. This method has 
a range of applicability deppending on the stellar rotation 
velocity (Suarez et al. 2005; Moya et al. 2010a). For in- 
clination angles around i = 36° (corresponding to a stel- 
lar rotation velocity of V mt ~ 60km s _1 (Kaye & Strass- 
meier 1998), the limit for the FRM to be accurately applied), 
the models fulfilling all the observational constraints have 
masses in two separate ranges of M=[1.32, 1.33], [1.44, 
1.45] Mq. The age of the system is constrained in two sep- 
arate ranges: [1123, 1625] Myr and [26, 430] Myr respec- 
tively (Fig. 2). A percentage of 16.7% of the models are 
in the range given in the literature. A consequence of this 
result is that, in the case of the youngest age range, the pre- 
dicted masses of the observed planets are [5,14] Mj up for 
the most luminous planets, and [3,13] Mj up for the less lu- 
minous one. The oldest age range, the most probable from 
the point of view of the present study, predicts masses for 
the three objects in the brown dwarfs domain (see Fig. 4 of 
Marois et al. 2008). 

The lack of an accurate determination of the inclination 
angle is the main source of uncertainty of the present study. 
This angle has not been unambiguously obtained up to now, 
and its value would say whether the results of this study are 
actually applicable, and then, it can provide a very accurate 
determination of the age and mass of the star. 

6 Conclusions 

The accurate determination of the mass, age and internal 
metallicity of a star hosting a planet or a planetary system is 
a necessary step to understand the formation and evolution 
of that planetary system. If the hosting star pulsates, the use 
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Table 2 Acceptable models depending on the physical constraints. In the complete procedure case, the first (second) 
mass range is linked with the first (second) age range shown. 



Constraint 


Mass 


Age 


% of models with [30, 160] 




in Mq 


in Myr 


in Myr 


HR position 


[1.25, 1.27], [1.32,1. 35], [1.40,1.48] 


[10, 2337] 


18.1 


Complete procedure 


[1.32, 1.33], [1.44, 1.45] 


[1123, 1625], [26,431 


)] 16.7 



of asteroseismology can provide these physical quantities 
with an accuracy hardly obtained with other techniques. 

The case of HR 8799 is a excellent example of this bene- 
fit provided by asteroseismology. In two recent works (Moya| 
et al. 2010a, b), the first comprehensive asteroseismologic 
study of the planetary system host HR 8799, a A Bootis star 
presenting 7 Doradus pulsations has been carried out. This 
asteroseismic work is specially important for the determina- 
tion of the internal abundances of this kind of stars, a previ- 
ous step to understand the physical mechanism responsible 
for the surface chemical peculiarities of the A Bootis group. 
On the other hand, an analysis of the age determination of 
the planetary system HR 8799 has also been done. The re- 
sults found in the literature are not conclusive, and the only 
valid argument to estimate the age of the star is that using 
its radial velocity and proper motion, but it is an estatistical 
argument needed of additional estimations. 

The selection of models fulfilling the spectroscopic ob- 
servations shows that there are no selected models with so- 
lar metallicity. This has an impact in the main assumption 
of the theory widely accepted to explain the A Bootis na- 
ture, i.e. that these stars have solar metallicity, whereas the 
observed abundances are due to surface phenomena. 

On the other hand, the main complementary argument 
to the kinematic age determination is the position of the star 
in the HR diagram. We have shown that, due to the A Boo- 
tis nature of HR 8799, that hides the internal metallicity of 
the star, this HR diagram position provides a range of ages 
[10,2337] Myr, a much wider than that estimated by other 
authors of [30,160] Myr (Marois et al. 2008). Only a small 
amount ( 1 8. 1 %) of models in our grid have ages in the range 
claimed in the literature. 

For inclination angles around i = 36°, the models ful- 
filling all the observational constraints have masses in two 
separate ranges of M=[1.32, 1.33], [1.44, 1.45] M , and 
the age of the system is constrained in two separate ranges: 
[1123, 1625] Myr and [26, 430] Myr respectively. A per- 
centage of 16.7% of the models are in the range given in the 
literature. This determination has an impact in the determi- 
nation of the mass of the objects observed orbiting around 
HR8799. 

A consequence of this study is the need for a precise 
determination of the inclination angle i, of the multicolour 
photometric amplitudes and phases of / 2 , and some infor- 
mation of m values through time-series if high resolution 
spectroscopy. These determinations would help to carry out 
a definitive selection of the models. In any case, the range 
of ages assigned to this star in the literature is unlikely to 



be the correct one. Only a stellar luminosity larger than that 
reported would allow young models with solar metallicity 
to fulfill all the observational constraints. 
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